AD-A270  935 


ON  PAGE 


*  •0*H  * 

D««>v  *< 


form  Appro^^ti 
OfytB  Afo  C704^cyss 


ir  '  •*Ow'  Of  ••xoc'\t  "t  *•'«  i  ’O'  ''9  c«<.l 

;f»*rn»o^  o*  *^a«»e'^9  O'  ♦''*  •vovr.  c*  ••'■k 

••♦•oc«*ne'*  Ve»*»<rk  *C'  Oor'»i«pf»>  fcppC'Tk  t;*5 

«ne  fc«.e3«'  fcpcwn'O^  ^»c>er  (C^C^'C  *H.  DC  ^CSC3 


1.  AGENCY  USE  only  (Lt^^t  Dunic; 


August  1993 


3  REPORT  TYPE  AND  DATES  COVERED 

final  report  Jan  1  1992  -  June  30  1 


A.  TITLE  AND  SUBTITLE 


Refractive  Index  Structure  in  the  Lower  Atmosphere 


6.  AUTHORIS) 

J  Wyngaard 
D  Thomson 


S.  FUNDING  NUMBERS 


DAAL03-92-G-0004 


93-24958 


73 


7.  PERFORMING  organization  NAME(S)  AND  ADDRESS(ES) 

Department  of  Meteorology 
The  Pennsylvania  State  University 
503  Walker  Bldg 
University  Park  PA  16802 


©Tte 

ELECT Ej 

0CC2  0.199* 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADORE 

U.  S.  Army  Research  Office 
P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


Ktri.^ni  nunnotn 


TO.  SPONSORING  /  MONITORING 
,  AGENCY  REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 

The  view,  opinions  and/or  findings  contained  in  this  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,-  policy,  or  decision,  unless  so  designated  bv  other  documentation. 


12*.  DISTRIBUTION /AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


12b.  DISTRIBUTION  CODE 


13.  abstract  (Maximum  200  words) 


I 

(  ■ 

€'■ 

( 


t 

C 


We  develop  a  local  scaling  theory  for  structure-function  parameters  C%  and  Cl  of 
refractive  index  and  velocity  fluctuations  in  the  stably  stratified  atmospheric  boimd- 
ary  layer.  The  theory  relates  them  to  the  local  T^alues  of  vertical  temperature  flux  w6 
and  shearing  stress  r.  The  theory  agrees  well  vdth  data  from  observations  and  from 
large-eddy  simulations.  The  effects  of  flow  unsteadiness,  baroclinity,  terrain  slope,  and 
breaking  giavity  waves  are  apt  to  cause  the  vertical  profiles  of  w6  and  r  to  vary  strongly 
from  case  to  case,  making  the  and  Cl  profiles  quite  variable  as  well,  but  according 
to  the  theory  the  stress  and  heat  flux  profiles  can  generally  be  diagnosed  from  those  of 
the  structure  parameters. 


U.  SUBJECT  TERMS 

refractive  index,  structure-function  parameters,  turbulence, 
stable  boundary  layer,  large-eddy  simulation 

15.  NUMBER  OF  PAGES 

33 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

IB.  security  CLASSIFICATION 

19.  security  classification 

OF  ABSTRACT 

20.  LIMITATION  OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

UL 

NSN  75AO-C-.-280-S500 


Stanoarc  Form  29E  (Rev  ;.S9) 

Cur  ANV’  StC 


REFRACTIVE  INDEX  STRUCTURE  IN  THE  LOWER  ATMOSPHERE 


Final  Report 


J.  Wyngaard,  D.  Thomson 
August,  1993 

U.  S.  Array  Research  Office 
Grant  Number  DAAL03-92-G-0004 


The  Pennsylvania  State  University 
Department  of  Meteorology 
University  Park,  PA 


Approved  for  Public  Release; 
Distribution  Unlimited 


Table  of  Contents 


A.  Statement  of  Problem  Studied . 3 

B.  Summary  of  Most  Important  Results . 3 

C.  Publications  and  Presentations . 6 

D.  Participating  Scientific  Personnel . 6 

E.  Bibliography . 6 

F.  Appendix . 8 


I  ioaasalon  for 

1  shs  GPAil 
1  DTIC  TAB 


Un'iniicuncod 


Justlf icatlo 

n 

- - - - — 

! 

i  ri  'i  o-.f.r  ;■  I'liti  on 
Av-.ilablllt 
jA’-cai  . 
Dl3t  ;  i'r'ia 

1^' 

/ 

y  Codes 
iuii/or 

I _  -  - -  - 

I 


2 


□  □ 


A.  Statement  of  Problem  Studied 


Theoretical  treatment  of  the  scattering  of  electromagnetic  and  acoustic  radiation 
usually  focuses  on  the  refractive  index  structure  function  (Tatarskii,  1971)  which  is  de¬ 
fined  by  _  _ 

[n(x,t)  -  n(x -l-r,<)]2  =  (n  -  n')2,  (1) 

where  n  is  the  fluctuating  refractive  index  at  point  x  and  n'  that  at  point  x  +  r.  The 
overbar  denotes  the  expected  value  (ensemble  average).  Many  wave- propagation  prob¬ 
lems  involve  this  structure  function  only  for  separations  r  =  |r|  in  the  inertial  range  of 
scales,  where  it  has  the  form 

(n  -  n')2  =  (2) 

We  call  Cy  the  refractive  index  stTuctuTe-function  varameteT. 

Fluctuations  in  the  refractive  index  for  microwave  radiation  are  caused  primarily  by 
fluctuations  in  temperature  and  absolute  humidity,  6  and  q,  respectively.  One  can  write 

n  =  aO  +  bq,  (3) 

where  the  coefficients  a  and  b  depend  on  the  wavelength  of  the  radiation,  Wjuigaard 
al.  (1978)  have  shown  that  this  implies  that  Cjf  may  be  expressed  in  terms  of  the 
structure-function  parameters  for  temperature  and  humidity,  Cj  and  Cq,  and  the  joint 
structure  peurameter,  Ctq- 

Cl  =  +  2abCTQ  +  (4) 

For  acoustic  radiation,  the  refractive  index  structure-function  parameter  depends  addi¬ 
tionally  on  Cl,  the  structure-function  parameter  for  the  velocity  field. 

The  purpose  of  our  research  was  to  illuminate  the  physics  of  the  maintenance  of 
these  structure-function  parameters  in  the  lower  atmosphere  and,  in  particular,  to  de¬ 
velop  a  theory  for  their  vertical  profiles  in  the  stably  stratified  boundary  layer. 

B.  Summary  of  Most  Important  Results 


1.  Local-scaling  theory  for  structure-function  parameters 


Using  a  combination  of  observational  and  large-eddy-simulation  data,  we  developed 
and  tested  a  similarity  theory  for  Cj.  and  Cl  in  the  stably  stratified  boundary  layer. 
This  theory  relates  these  structure-function  parameters  to  the  vertical  temperature  flux 
and  the  shearing  stress,  basic  structural  variables  in  the  stable  boundary  layer.  Our  the¬ 
ory  uses  the  “local  scaling”  hypothesis  for  stably  stratified  shear  flows.  In  this  hypothe¬ 
sis  the  length,  velocity,  and  temperature  scales  for  the  turbulence  are  taken  as 


kgw6  ’ 


(5) 


where  t{z)  and  wB{z)  are  the  local  values  of  kinematic  stress  and  temperature  flux,  re¬ 
spectively.  The  basic  hypothesis  is  that  dependent  variables  made  dimensionless  with 
these  sczdes  are  universal  functions  of  zf k.  In  effect,  this  generalizes  Monin-Obukhov 
similarity  to  the  entire  stable  boundary  layer  by  replacing  the  the  Monin-Obukhov 
length  scale  Tj  with  A. 

The  local  scaling  hypothesis  has  roots  both  in  the  behavior  observed  in  the  1968 
Kanseis  surface-layer  experiments  and  other  field  programs  (Wyngaard,  1973)  as  well  as 
in  Nieuwstadt’s  (1984)  studies  of  a  set  of  second-moment  equations  for  the  stable  bound- 
axy  layer.  The  Kansas  data  for  wind  and  temperature  gradients,  velocity  and  temper¬ 
ature  variances,  and  dissipation  rates  at  large  zJL  (i.e.,  imder  very  stable  conditions  in 
the  surface  layer)  are  consistent  with  local  scaling.  Nieuwstadt  found  that  velocity  and 
temperature  statistics  from  the  200-m  tower  at  Cabauw  were  also  in  good  agreement 
with  local-scaling  predictions. 

The  local-scaling  prediction  for  is 


T} 


(6) 


where  /  is  a  universal  function  to  be  determined.  In  the  very  stable  surface  layer  the 
Kansas  data  yield  the  M-0  asymptote 


T2 


20, 


(7) 


(Wynganrd  et  al.,  1971).  Thus,  if  local  scaling  exists  we  expect  the  made  dimension¬ 
less  with  local  scales  to  be  constant  at  a  magnitude  of  about  20  in  the  outer  regions  of 
the  stable  boundary  layer  as  well. 

We  tested  this  prediction  with  numerical  data  from  the  Meison  and  Derbyshire 
(1990)  large-eddy  simulation  (LES)  of  the  stable  boundary  layer  and  with  observational 
data  from  the  1973  Minnesota  experiments  (Caughey  et  al.,  1979).  The  LES  domain  wets 
500  m  long  in  the  direction  of  geostrophic  wind,  300  m  wide,  and  1000  m  deep.  This  do¬ 
main  had  40  X  32  X  62  grid  points.  For  analysis  of  these  numerical  data  we  used  the 
Cray-YMP  computer  facility  at  the  National  Center  for  Atmospheric  Research  (NCAR). 
In  all  cases  we  used  spatial  averages  over  horizonted  planes  instead  of  ensemble  averages. 

We  fotmd  that  over  most  of  the  stably  stratified  layer  the  dimensionless  from 
the  Minnesota  measurements  is  in  the  range  10-20.  The  dimensionless  Cf  from  the 
Mason- Derbyshire  LES  results  is  in  the  range  12-16  over  the  stable  PBL.  These  show 
fairly  good  agreement  with  the  surface-layer  eisymptote  of  20  found  in  the  Kansas  data. 
For  the  dimensionless  C^,  the  Minnesota  data  are  in  the  range  10-15  while  the  LES  data 
yield  a  value  of  about  7.  The  Kansas  surface-layer  asymptote  is  about  10,  in  fairly  good 
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agreement  with  these  other  results.  Thus,  our  results  support  the  notion  of  local  scaling 
of  structure-function  parameters  in  the  stably  stratified  boundary  layer. 

Our  results  also  show  that  the  key  governing  parameters  in  local  scaling — the  tem¬ 
perature  flux  w6{z)  and  shearing  stress  t{z) — are  very  sensitive  to  imsteadiness  of  the 
boundary-layer  flow,  baroclinity,  terrain  slope,  and  breaking  gravity  waves.  As  a  result, 
the  vertical  profiles  of  w9  and  r  are  apt  to  vary  strongly  from  case  to  case,  depending 
on  the  magnitudes  of  these  effects.  It  follows  from  locJil  scaling  that  the  profiles  of  Cj- 
and  variable  as  well.  The  Cj-  profile  is  particularly  variable;  it  can  either  increase  or 
decrejise  with  height,  as  found  experimentally  by  Cuijpers  and  Koshiek  (1989). 

The  appendix  contains  a  full  report  on  these  results. 

2.  Structure-function  parameters  in  the  convective  boundary  layer 

We  also  began  a  study  of  the  structure-function  parameters  for  temperature  and 
humidity  in  the  convective  boundary  layer.  This  study  was  motivated  in  part  by  a  re¬ 
quest  from  the  Microwave  Remote  Sensing  Laboratory  (MIRSL)  at  the  University  of 
Massachusetts  for  information  on  €%  fields  in  the  boundary  layer.  They  need  such  data 
for  design  studies  of  their  Turbulent  Eddy  Profiler  (TEP)  radar  being  developed  imder  a 
URI  grant  from  ARO. 

As  an  ensemble-average  quantity,  the  structure-fimction  parameter  does  not  fluctu¬ 
ate.  In  propagation  applications,  however,  it  can  be  necessary  to  define  a  “local”  struc¬ 
ture  function  parameter  which  does  have  a  random  component.  The  traditional  temper¬ 
ature  structure-function  parameter  is  proportionzd  to  the  inertial-subrange  level  of  the 
temperature  spectrum,  and  so  by  the  similarity  theory  for  that  spectrum  can  be  written 
(Wyngaard  et  al.,  1971) 

(8) 

Here  e  and  e$  are  the  (ensemble-average)  molecular  destruction  rates  of  turbulent  kinetic 
energy  and  temperature  variance,  respectively,  and  /?  is  the  spectral  constant.  A  local 
temperature  structure-function  parameter  is 

=  (9) 

where  the  tilde  denotes  quantities  averaged  over  a  loced  spatial  volume  rather  than  an 
ensemble  of  realizations. 

With  L.  Joel  Peltier,  who  is  supported  by  our  URI  grant  from  ARO,  we  are  us¬ 
ing  LES  data  kindly  provided  by  Chin-Hoh  Moeng  of  NCAR  to  analyze  the  behavior 
of  the  local  structure-function  parameter  fields  for  temperature  and  humidity  (defined 
by  averaging  over  the  LES  grid  volume)  in  the  convective  botmdary  layer.  The  results 
show  that  the  vertical  profiles  of  their  area  averages  (which  should  correspond  to  the 
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traditional,  ensemble-average  parameters)  agree  well  with  aircraft  data  (Wyngaard  and 
LeMone,  1980).  However,  the  local  parameters  fluctuate  considerably;  one  measure  of 
their  variability  is  the  ratio  of  variance  to  squared  mean. 


A'  = 


(Cl? 


(10) 


where  the  overbar  denotes  averaging  over  the  horizontal  plane.  The  LES  data  indicate 
that  K  is  of  the  order  of  2-3  in  the  convective  boundary  layer.  We  have  begun  an  ex¬ 
tended  series  of  acoustic-sounder-based  Cj  measurements  at  our  Rock  Springs,  PA  field 
site,  from  which  we  will  determine  K  experimentally. 

This  work  is  continuing  under  our  URI  grant  from  ARO. 


C.  Publications  and  Presentations 

We  have  no  formal  publications  on  this  work  during  the  period  of  the  grant,  but  the 
appendix  is  the  draft  of  a  manuscript  “Similarity  of  structure-fimction  parameters  in  the 
stably  stratified  boundary  layer”  that  will  be  submitted  for  publication  shortly. 

Dermis  Thomson  presented  an  invited  paper  “Wind-field  and  turbulence  sensing”  at 
the  IEEE  Topical  Symposium  on  Combined  Optical-Microwave  Earth  and  Atmosphere 
Sensing,  22-25  March,  1993  in  Albuquerque.  He  then  traveled  to  White  Sands  Missile 
Range  to  spend  Meirch  26  briefing  personnel  in  the  Battlefield  Environment  Directorate 
(BED)  on  current  research  in  boundary-layer  meteorology  that  relates  to  tactical  envi¬ 
ronmental  support  systems.  A  number  of  opportimities  for  continuing  cooperative  re¬ 
search  including  FM-CW  radar  and  atmospheric  acoustical  propagation  studies  were  also 
discussed.  Doug  Brown  of  BED  organized  and  hosted  Prof.  Thomson’s  visit. 


D.  Participating  Scientific  Personnel 

Faculty:  Dennis  Thomson,  John  Wyngaard. 

Graduate  students;  Robert  Edsall,  Jeffrey  Hare,  Branko  Kosovic. 
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APPENDIX 


SIMILARITY  OF 

STRUCTURE-FUNCTION  PARAMETERS 
IN  THE  STABLY  STRATIFIED  BOUNDARY  LAYER 

J.  Wyngaard  and  B.  Kosovic 
Department  of  Meteorology 
Penn  State  University 


Abstract 

Published  data  on  the  structure-function  parameters  and  of  temperature  and 
velocity,  respectively,  follow  Monin-Obukhov  (M-0)  simileirity  in  the  stable  surface  layer. 
In  the  outer  layer,  Cf  and  Cl  data  from  the  1973  Minnesota  experiments  and  from  the 
Mason-Derbyshire  large-eddy  simulation  study  support  Nieuwstadt’s  local- scaling  hy¬ 
pothesis.  The  data  suggest  a  smooth  transition  from  M-0  similarity  to  local  similarity. 

The  similzuity  scales  in  the  outer  layer  depend  on  the  local  values  of  turbulent  stress 
and  buoysincy  flux,  suggesting  that  remote-sensing  measurements  of  and  Cl  could  be 
used  to  infer  profiles  of  stream  and  buoyancy  fiux.  The  lo  al- scaling  hypothesis  implies, 
through  the  mean-flow  conservation  equations,  that  the  vertical  profiles  of  structure- 
function  parameters  eu’e  sensitive  to  unsteadiness,  baroclinity,  terrain  slope,  and  breaking 
gravity  waves.  We  show  that  as  a  result  the  profile,  in  particular,  can  vary  consider¬ 
ably  from  case  to  case.  This  is  observed. 


8 


1.  Introduction 


Many  boundary-layer  researchers  will  recall  the  impact  of  the  McAllister  et  al. 

(1969)  paper.  They  had  pulsed  an  array  of  loudspeakers  and  detected  the  backscattered 
energy  with  the  stime  array.  By  cleverly  recording  the  signal  on  a  facsimile  receiver,  they 
produced  striking  time-height  displays  showing  how  convective  eddies  rising  from  the 
heated  surface  in  the  early  morning  break  up  the  nocturnal  boimdary  layer.  Use  of  the 
acoustic  echo  sounder,  as  it  was  soon  called,  spread  rapidly;  before  long,  researchers  were 
studying  facsimile  returns  from  shallow  but  distinctly  turbulent  nocturnal  PBLs. 

At  that  time  the  expected  height  h  of  the  nocturnal  boundary'  layer  was  a  subject 
of  speculation  and  controversy  (DeardorfF,  1972).  Similarity  theories  (e.g.,  Clarke,  1970) 
held  that  h  satisfied  the  equilibrium  relation 

“(T)'(7i)’ 

where  c  ~  0.3,  u,  is  the  friction  velocity  (the  square  root  of  the  kinematic  surface  stress), 
/  is  the  Coriolis  parameter,  and  L  is  the  Monin-Obukhov  length;  the  fimction  F  is  unity 
at  neutral  and  decreases  under  stable  stratification.  It  was  not  clear,  however,  whether 
the  nocturnal  boundary  layer  could  reach  this  equilibrium.  The  behavior  of  conventional 
soundings  of  temperature  and  wind  had  led  Blackadar  (1957)  to  conclude  that  h  typi¬ 
cally  increases  during  the  night  from  its  early-evening  minimum.  DeardorfF  suggested 
that  the  acoustic  echo  soimder  would  be  able  to  provide  direct  data  on  h{t)  in  order  to 
study  this  question. 

DezirdorfF  concluded  on  the  basis  of  his  one-dimensional  model  calculations  that  h  is 
generally  time  dependent  and  cannot  be  diagnosed  from  a  steady-state  similarity  relation 
such  as  (1).  [He  had  just  cau'ried  out  time-dependent,  fine-mesh,  three-dimensional  cal¬ 
culations  of  the  neutral  and  convective  PBLs;  his  attempts  to  extend  these  to  the  stable 
case  failed  because  most  of  the  turbulence  became  subgrid-scale.  This  approach,  which 
we  now  call  Izu'ge-eddy  simulation,  or  LES,  was  first  applied  successfully  to  the  stable 
case  by  Mason  and  Derbyshire  (1990).] 

Rather  than  prescribing  eddy-difFusivity  profiles,  as  DeardorfF  had  done,  Delage 
(1974)  used  eddy  difFusivities  based  on  predicted  turbulent  kinetic  energy.  His  one¬ 
dimensional  model  of  the  nocturnal  boundary  layer  reached  a  steady  state  after  several 
hours  with  a  realistic  rate  of  surface  cooling.  Wyngaard  (1975),  integrating  a  full  set  of 
second-moment  equations,  had  similar  findings. 


Brost  and  Wyngaard  (1978)  used  a  simplified  form  of  Wyngaard's  second-moment 
model  to  calculate  the  structure  and  evolution  of  the  stable  boundary  layer,  \^’ith  a  con¬ 
stant  surface- cooling  rate  it  .  ^uld  reach  an  essentially  steady  state  within  a  few  hours. 
Nieuwstadt  (1984)  shew  that  the  Brost-Wyngaard  equations  have  steady  solutions 
that  display  a  “local”  similarity  that  is  an  extension  of  Monin-Obukhov  (M-0)  simi¬ 
larity.  In  M-0  similarity  the  friction  velocity  a,  and  the  surface  temperature  flux  Bu't 
define  the  Obukhov  length  L: 

L  = - (2) 

H9lT)Bw, 

Here  k  is  the  von  Karman  constant  and  g/T  is  the  buoyancy  parameter.  The  counter¬ 
part  of  L  in  local  similarity  is  the  local  Obukhov  length  A 

3/2 

A  - - =,  (3) 

k{gfT)Bw 


where  t{z)  is  the  magnitude  of  the  local  kinematic  stress  (i.e.,  =  uw^  -f  vw^) 

and  Bw{z)  is  the  local  temperature  flux.  Under  Nieuwstadt’s  “loceil  scaling”  hypothe¬ 
sis,  quantities  made  dimensionless  with  r,  Bw,  and  A  depend  oiJy  on  z/A.  Nieuwstadt 
showed  that  obser\'ations  taken  in  the  quasi-steady  nocturnal  PBL  at  Cabauw,  The 
Netherlands,  support  this  hypothesis. 

Nieuwstadt  pointed  out  that  while  local  scaling  may  be  a  theoretically  satisfac¬ 
tory  description  of  stable  turbulence,  it  is  not  very  well  suited  for  practical  applications, 
which  usually  require  verticeJ  profiles  of  turbulence;  local  scaling  does  not  give  profiles. 
However,  one  can  obtain  profiles  by  using  local  scaling  as  a  closure  in  model  equations 
for  the  stable  boundzury  layer.  One  finds  that  in  generzil  the  turbulence  profiles  depend 
on  the  time  history  of  the  boundary  layer.  As  Wyngziard  (1988)  states,  "...  a  simple, 
diagnostic  similzirity  framework  for  the  structure  of  the  nocturnal  boundary  layer  seems 
not  to  exist.  What  we  have  described  instead  is  a  relatively  simple  model  of  its  dynam¬ 
ics  _  (Its)  solution  . . .  gives  the  structure,  which  then  depends  on  the  historj'  of  the 

flow  and  boundary  conditions.” 


2.  The  Structure-function  P'xrameters 
A.  Refractive  index 

Most  theoretical  analyses  of  the  scattering  of  electromagnetic,  acoustic,  and  opti¬ 
cal  radiation  by  turbulence  use  the  refractive  index  structure  function  (Tatarskii,  1971) 
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which  is  defined  by 

[n(x,i)  -  n(x  +  r,#)]^  =  {n  -  n')^,  (4) 

where  n  is  the  fluctuating  refractive  index  at  point  x  and  n'  that  at  point  x  +  r.  The 
overbar  denotes  the  ensemble  average.  Many  problems  of  practical  importance  involve 
this  structure  function  only  for  separations  r  =  |rj  in  the  inertial  range  of  scales,  where  it 
has  the  form 

{n-n')^  =C%r'^^^.  (5) 

Fluctuations  in  the  refractive  index  for  electromagnetic  radiation  are  caused  pri¬ 
marily  by  fluctuations  in  temperature  and  absolute  humidity,  B  and  q,  respectively.  For 
fluctuation  levels  t3'pical  of  the  atmospheric  boimdary  laj’er  one  can  write 

n  =  aB  +  bg,  (6) 

where  the  coefficients  a  cmd  b  depend  on  the  wavelength  of  the  radiation.  Wjmgaard 
et  al.  (1978)  have  shown  that  this  implies  that  may  be  expressed  in  terms  of  the 
structure-function  parameters  for  temperature  and  humidity,  and  Cq,  and  the  joint 
parameter,  Ctq- 

Cl,  =  +  2abCTQ  +  b^C^Q.  ( 7) 

(For  acoustic  radiation  Cl,  depends  also  on  the  turbulent  velocity  field.)  Functioned 
forms  for  the  coefficients  in  (7)  are  summarized  by  Burk  (1979)  for  acoustic,  microwave, 
eind  opticcd  radiation.  Gossard  (1988)  heis  discussed  the  maintenance  of  the  joint  param¬ 
eter  Ctq  and  its  role  in  radar  scattering. 

These  structure-fxmction  parameters  are  directly  proportional  to  the  spectra  of  tem¬ 
perature  and  humidity,  and  their  cospectrum,  respectively,  in  the  inertial  subrange.  Sim¬ 
ilarity  arguments  due  to  Kolmogorov,  Obukhov,  and  Corrsin  (see  Teimekes  zmd  Lum- 
ley,  1972)  indicate  that  these  inertial-range  spectral  levels  depend  only  on  the  molecular 
destruction  rates  e,  Xff,  Xq,  and  X0q  of  turbulent  kinetic  energy,  potential  temperature 
and  humidity  A'ariances,  zmd  their  covariance,  respectively'.  In  the  inertial  range  the  one- 
dimensioned  spectra  behave  as  (Wyngaard  and  LeMone,  1980): 

$r(«i)  =  2(27r)-Jr(5/3)sin(7r/3)C^Kr"/" 

=  0.25C7-Kf®'^^  = 

Cotq{ki)  =  ^.2bCTQftx^^^  = 

$Q(/ca)  =  0.25CjKr"/"  = 


(8) 


where  and  71  are  the  one-dimensional  spectral  constants.  Our  normalization  in  (8) 
is  such  that  the  integral  over  the  half-line  is  the  co\-ariance.  Following  Wyngaard  and 
Lemone  (1980),  we  take  =  7j  =  0.4  and  recast  (8)  as 

Ctq  =  (9) 

Cl  = 

B.  Velocity 

As  discussed  by  Kaimal  (1973),  the  inertial-subrange  form  of  the  velocity  structure 
function  for  uj,  for  example,  is 

[ui(x,t)  -  ui(x  +  r,0]^  =  (10) 

where  is  called  the  velocity  structure-function  pareuneter.  Kaimal  presents  the  spec¬ 
tral  relations  analogous  to  (8).  For  our  purposes  we  need  only  the  result 

=  4aic2/3  ~  2f2/^  (11) 

where  o  j  ~  0.5  is  the  one-dimensional  spectral  constant  for  velocity.  When  acoustic 
radiation  is  scattered  at  angles  other  than  180  deg,  can  be  the  dominant  factor  in  the 
scattering  cross  section. 

3.  Similarity  Structure 

Eqs.  (9)  and  (11)  tie  the  structure-function  parameters  to  the  rates  of  molecular  de¬ 
struction  of  the  temperature-humidity  covariances  and  of  turbulent  kinetic  energ>'.  We 
can  relate  these  destruction  rates  to  the  dynamics  of  the  turbulence  through  the  second- 
moment  equations.  We  derive  these  formally  by  using  the  mean-fluctuating  decomposi¬ 
tion  for  cdl  random  fields  (Tennekes  and  Lumley,  1972).  We  denote  these  fields  with  a 
tilde,  using  upper  and  lower  case  symbols  for  their  mean  and  fluctuating  parts,  respec¬ 
tively: 

Ui  =  Ui-\-Ui-,  p  =  P  +  p;  f  =  Q  +  6-,  q  =  Q  +  q.  (12) 

We  will  use  the  mixed  notation  Uj  =  (uj,  U25  ^3)  =  (u,  v,  w);  Xi  =  (xi,  r2,  ^3)  =  (2:,  y,  z). 
Let  us  assume  that  conditions  are  locally  homogeneous  in  the  horizontal,  so  that  mean 
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advection  is  negligible  in  the  second-moment  equations.  The  potential  temperature  vziri- 
ance  budget  is  (Wyngaard  and  Cote,  1971): 


d9^ 

dt 


az 


+ 


d$^w 

dz 


X 


(13) 


The  humidity  variance  and  temperature-humidity  covariance  budgets  are  (Wyngaard  et 
al.,  1978) 


where  p  is  the  mean  density.  FinaJly,  the  turbulence  kinetic  energy*  (TKE)  equation  is 


d  Ui'Ui 

'm~2~ 


_ dU  __dV 

uw-z - 1-  VW-Z - h 

az  az 


d  WUiUi 

Tz  2 


1 

p  dz 


+  ^9w  -  e, 


(16) 


where  repeated  indices  are  summed.  We  are  assuming  here  that  the  humidity  is  suffi¬ 
ciently  small  that  it  does  not  contribute  significantly  to  buoyancy;  if  this  is  not  the  case, 
one  can  use  the  same  equations  but  interpret  T  as  the  virtual  potential  temperature 
(Lximley  and  Panofsky,  1964), 

It  seems  generally  agreed  (Panofsky  and  Dutton,  1984)  that  to  a  good  approxima¬ 
tion  in  the  quasi-steady,  locally  homogeneous  surface  layer  the  budgets  of  scalar  variance 
and  turbulent  kinetic  energy'  reduce  to  a  local  balance  between  production  and  molec¬ 
ular  destruction.  The  9q  budget  has  not  received  as  much  attention,  but  Wyngziard  et 
al.  (1978)  foimd  that  in  the  unstable  surface  layer  it  does  have  the  same  local  balance. 
Since  our  experience  with  these  second-moment  budgets  generciUy  indicates  that  the  tur¬ 
bulent  transport  (i.e.,  the  third-moment  divergence)  terms  are  least  important  under  sta¬ 
ble  conditions,  it  is  plausible  that  this  production-destruction  balance  for  the  9q  budget 
exists  under  stable  conditions  as  well.  Thus,  we  write  (13)-(15)  as 


Xe 
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(17) 


The  TKE  equation  (16)  becomes 


_ dU  dV  g  s- 

t  =  —UW— - VXV-—-  +  ~tfw. 

dz  dz  To 


(18) 


A.  Surface-layer  similarity 

We  define  surface-layer  temperature  and  humidity  scales  T*  and  9,, 


T,  =  - 


Bwt 


9*  = 


u. 


(19) 


where  is  the  surface  humidity  flux.  Under  the  M-0  hypothesis  in  the  quasi-steady, 
locally  homogeneous  surface  layer  mean  gradients  made  dimensionless  with  these  scales 
and  height  2 


kz  dQ 

~  fill' 

_  ^  _a  fQ\ 
q*^dz\p)' 


(20) 


are  universal  functions  of  z/L.  It  is  not  clear  whether  (f>h  and  4>g  differ;  we  will  adopt  the 
suggestion  of  Panofsky  and  Dutton  (1984)  that  they  be  assumed  equal.  In  the  surface 
layer  the  fluxes  of  temperature,  humidity,  and  momentum  are  essentially  equal  to  their 
surface  values,  so  under  the  M-0  hypothesis  we  can  write  (17)  as 


2  2  2  u,T»q» 

—  U  0^1  X?  t  7  0/t' 


(21) 


We  choose  the  z-axis  along  the  mean  wind,  so  that  Ui  =  [1/(2),  0,0]  and  the  M-0 
function  for  mean  wind  shear  is 

i-r  riTT 

(22) 


,  kzdU 

'm  -  ej 

U*  dz 


We  can  write  the  dissipation- rate  expression  (18)  as 


Combining  (9),  (21),  and  (23)  then  yields  M-0  expressions  for  the  refractive  index 
structure-function  parameters: 

2  3.2  Tl  ^  1/3 

—  ^2/3  ^2/3*^*  ’ 


(23) 
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(24) 


That  Tor  velocity  is 


3.2  ql 


-1/3 


Ch  = 


~  ^2/3 


(25) 


There  are  two  important  aspects  of  the  expressions  in  (24).  First,  they  predict  that 
the  structure-function  parameters  made  dimensionless  with  the  surface-layer  scales  vary 
identically  with  the  stability  p<irameter  z/L.  Second,  the  structure  parameter  Ctq  is 
not  positive  definite,  as  are  the  other  two;  its  sign  is  that  of  T»q^ — i.e.,  that  of  the  prod¬ 
uct  of  the  surface  fluxes  of  temperature  and  humidity.  Under  stably  stratified  conditions 
the  surface  temperature  flux  is  negative;  we  would  normally  expect  the  surface  humidity 
flux  to  be  positive,  mahing  Ctq  negative. 

Of  the  three  refractive-index  structure-function  parameters  only  Ct  has  been  mea¬ 
sured  extensively.  The  Kansas  results  (Businger  et  al.,  1971)  indicate  for  stable  condi¬ 
tions 

<f>h  =  0.74 +  A.7 z/L-,  <i>m  =  1.0  4- 4.7 z/L,  (26) 

so  that  with  the  Kansas  value  of  0.35  for  k  the  M-0  form  (24)  for  Ct  becomes 

_  4.8(1 +6.35Z/L) 

t;  (1  +  z.7z/Ly/^ '  '  •’ 

Figure  1  shows  Cj-  data  from  the  1968  Kansas  experiments  (Wyngaard  et  al.,  1971),  the 
1973  Minnesota  experiments  (Caughey  et  al.,  1979),  8ind  from  experiments  reported  by 
Foken  and  Kretschmer  (1990).  Eq.  (27)  fits  those  results  fairly  well. 


B.  The  stable  surface-layer  asymptote 

Under  very  stable  conditions  (large  z/L)  our  expression  (27)  for  Ct  becomes 

<^2  r2/3 

(28) 

Eq.  (25)  for  C^  becomes,  using  (26), 


u 


2 


~  9.0{z/Lfl\ 


(29) 
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These  can  be  rewritten  as 


Cf.20^, 


Cl  ~  9.6 


12/3- 


(30) 


A  “local  less  stratification”  interpretation  of  (30)  (Wyngaard,  1973)  is  that  under 
very  stable  conditions  turbulent  eddy  sizes  are  restricted  by  stability  and  z  loses  signif¬ 
icance  as  a  length  scale  for  the  turbulence.  Then  the  onlj'  length  scale  Z,  the  only  tem¬ 
perature  scale  is  T,,  and  the  only  velocity  scale  is  u,,  so  that  and  must  scale  with 
and  respectively.  Both  the  flux  and  gradient  Richardson  numbers  ap¬ 

proach  constants  interpreted  as  their  limiting  values  for  stable  stratification. 


C.  Outer- layer  similarity 

Let  us  briefly  review  the  physics  underlying  the  notion  of  local  scaling  of  turbulence 
in  the  stably  stratified  boundary  layer  and  the  evidence  supporting  it. 

Progress  in  turbulence  has  traditionally  been  paced  by  observational  studies  (Brad¬ 
shaw,  1972).  Because  observ'ations  are  usueJly  more  difl5cult  and  expensive  in  geophys¬ 
ical  flows  than  in  the  laboratory,  our  understanding  of  the.  atmospheric  boundcu-y  layer 
has  lagged  that  of  engineering  boundary  layers.  When  the  requisite  micrometeorologi- 
cal  turbulence  sensors  and  data-acquisition  systenas  became  available  by  the  1960s,  re¬ 
searchers  studied  in  detail  the  turbulence  structure  of  the  surface  layer  and  the  transport 
processes  within  it.  The  1968  Kansas  experiments  (Haugen  et  al.,  1971),  for  example, 
documented  this  structure  over  a  wide  range  of  stability  conditions.  Such  experiments 
at  ideal  sites  showed  that  under  stable  conditions  the  second-moment  budgets  displayed 
a  characteristic  state  of  local  equilibrium  in  which  the  time-change  and  turbulent  trans¬ 
port  (third-moment  flux  divergence)  terms  were  small  and  the  production  and  destruc¬ 
tion  rates  were  in  balance,  as  exhibited  in  Eqs.  (17)  and  (18).  Furthermore,  the  data  be¬ 
haved  as  if  L  rather  than  z  determines  the  length  scale  of  the  turbulence.  These  findings 
stood  in  striking  contest  with  those  for  the  unstable  surface  layer. 

Numerical  models  of  the  stably  stratified  boundary  layer,  using  closures  rooted  in 
this  interpretation  of  the  surface-layer  physics,  gradually  appeared  (Delage,  1974;  Brost 
and  Wyngaard,  1978;  Nieuwstadt,  1984).  They  might  have  seemed  bold  extrapolations, 
but  there  were  indications  that  they  could  accoimt  fairly  well  for  behavior  observed  in 
at  least  some  experiments  (Caughey  et  al.,  1979;  Nieuwstadt,  1984).  However,  others 
showed  features  such  as  flow  intermittency  under  very  stable  conditions  (Kondo  et  al., 
1978)  and  gravity-wave  effects  (de  Baas  and  Driedonks,  1985)  that  were  not  accoimted 
for  or  reproduced  by  munerical  models. 


When  the  stably  stratified  boundary  layer  yielded  to  large-eddy  simulation  (Ma¬ 
son  and  Derbyshire,  1990),  a  new  avenue  to  studying  its  structure  and  dynamics  was 
opened.  Derbyshire  (1990)  soon  suggested  that  a  key  parameter  is  BojG^f,  where  Bq  is 
the  surface  buoyancy  flux  and  G  is  the  magnitude  of  the  geostrophic  wind.  He  argued 
that  when  this  parameter  is  less  than  a  threshold  vaJue  the  layer  can  exist  in  a  quasi¬ 
steady  state.  Mzison  and  Derbyshire  stated  that  “The  potentially  nonlocal  character  of 
some  wave  motions  does  not  seem  to  be  strong  enough  to  invalidate  local  scaling  .... 
Arguments  given  by  Derbyshire  (1990)  show  that  this  is  to  be  expected  if  the  waves  are 
generated  purely  by  turbulence,  though  topographic  waves  may  of  course  behave  differ¬ 
ently.”  In  the  quasi-steady  state  and  without  topography,  the  Mason-Derbyshire  results 
indicate  broad  agreement  with  the  local-scaling  hypothesis  for  the  nocturnal  boundary 
layer. 

Under  the  local  scaling  hypothesis  we  cein  write  for  the  outer  layer 


Cr  =  ^  /.(^/A),  Cl  =  ^  Mz/A).  (31) 

where  /j  and  /2  are  functions  to  be  determined.  The  vanishing  significance  of  z  under 
strong  stratification  implies  that  they  approach  constants  in  the  limit  2/ A  -4  00.  The 
flux  and  gradient  Richardson  numbers  are  also  constant  at  their  limiting  values.  The 
Brost-Wyngaard  model  gives  in  this  stable  limit 


/i  = 


6.2u)2 
£^5/6  . 


1I/6 


f2  = 


0.549 


.=2 


’  Kl 


/3 


(32) 


where  Kh  is  the  eddy  diffusivity  for  temperature,  is  twice  the  turbulent  kinetic  energy 
per  unit  mass,  and  a  tilde  denotes  nondimensionalization  with  the  local  scales  dw,T^  and 
A.  Each  of  these  dimensionless  dependent  variables  is  a  function  of  2/A.  Nieuwstadt's 
numerical  solutions  show  that  as  2/ A  increases  they  approach  constant  values  of 


Kh  ^  C.08,  q  ~  2.9,  ~  i.7, 


(33) 


so  that  from  (32) 


/i  ~  19,  /2  ^  9, 


and 


tA2/3’ 


(34) 

(35) 
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These  agree  well  with  the  surface-layer  asymptotes  in  Eq.  (30),  because  the  closure 
constants  in  the  Brost-Wyngaard  model  were  chosen  so  that  its  results  were  compati¬ 
ble  with  the  stable  limit.  This  is  consistent  with  Sorbjan’s  (1986)  hypothesis  that  for  a 
given  variable  the  forms  of  the  M-0  and  outer-layer  similarity  functions  are  the  seime. 

We  tested  the  predictions  (35)  of  the  local-scaling  hypothesis  with  LES  results  and 
with  experimental  data.  The  results  from  the  Mason-Derbyshire  (1990)  LES  study  of  the 
stably  stratified  boundary  layer  indicate  that  made  dimensionless  with  local  scales 
approaches  an  asymptote  of  about  12  (Figure  2),  and  that  so  nondimensionalized  has 
an  asymptote  of  about  7  (Figure  3).  These  asymptotes  agree  fairly  well  with  those  of 
(35). 

The  experimental  data  were  talcen  at  0.3  <  z/h  <  0.75  in  the  1973  Minnesota  ex¬ 
periment  (Caughey  et  al.,  1979).  Unfortunately  only  uw  data  are  available;  the  lateral 
stress  was  contaminated  by  balloon  motion.  Thus,  in  calculating  the  local  scaling  pa¬ 
rameter  r  we  neglected  the  contribution  of  vw,  which  is  difficult  to  justify  near  h.  The 
Minnesota  Cj-  and  data  so  scaled  are  plotted  in  Figures  4  zmd  5,  respectively.  In  gen¬ 
eral  they  agree  well  with  the  local-scaling  results  (35)  and  the  M2ison-Derbyshire  LES 
results. 


4.  Vertical  Profiles 


From  (30)  and  (35)  we  can  write  for  most  of  the  stable  bovmdary  layer — from  the 
outer  edge  of  the  surface  layer  through  the  outer  layer — 


Cliz)  ~  20 


\  t)  r2  ’ 


cUz)~9[kiy'‘\ew\^i\ 


(36) 


Eq.  (36)  shows  that  the  profile  depends  on  the  profiles  of  kinematic  stress  mag¬ 
nitude  t(z)  and  temperature  flux  dw(zy,  we  will  see  that  the  balance  here  is  so  delicate 
that  Cj’  can  either  increase  or  decrease  with  height.  By  contrast,  the  profile  of  de¬ 
pends  only  on  the  8w  profile,  so  that  we  expect  it  to  decrease  monotonically  with  height. 
Let  us  now  consider  the  behavior  of  the  r  and  8w  profiles  in  more  detail. 


A.  Steady  conditions 

In  the  simplest  steady  case  with  horizontal  homogeneity  and  r-independent 
geostrophic  wind  components  Ug  and  Vg,  the  equations  for  the  vertical  gradients  of  mean 
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wind  and  temperature  reduce  to 


d^uw  dV  d^vw 
5^2  ~d^ 


(37) 

(38) 


Using  the  boimdary  conditions  that  wB  =  6wg  at  2  =  0  and  0w  =  0  a.t  z  =  h,  (38)  yields 
the  steady-state  6w  profile 

6u'  =  ^u’4(l  —  z/h).  (39) 


Nieuwstadt  (1984)  showed  also  that  a  stress  closure  based  on  local  scaling  gives  a  solu¬ 
tion  of  (37)  whose  r  profile  is 

T  =  ul{l  —  z/h)^^^.  (40) 

It  follows  from  (36)  that  the  steady  profiles  of  Cj>  and  from  the  outer  edge  of  the 
surface  layer  through  the  outer  layer  are 

Ct  ^  cl  ~  9^(1  -  z/hf/K  (41) 


Eq.  (41)  predicts  that  Cj-  increases  with  height,  but  the  observations  of  Caughey  et 
al.  (1979)  and  Cuijpers  and  Koshiek  (1989)  show  that  it  decreases.  The  prediction  of  an 
increasing  profile  results  from  the  steady-state  balance  between  the  competing  effects 
of  dw,  which  causes  to  decrease  with  height,  and  r,  which  causes  it  to  increase.  Let 
us  now  exEunine  other  factors  that  can  influence  the  shapes  of  the  9w  and  r  profiles  and, 
hence,  change  this  balance. 

B.  Influences  on  the  structure-parameter  profiles 


The  stress  and  temperature-flux  profiles  are  coupled  to  the  mean  horizontal  mo¬ 
mentum  and  mean  potential  temperature  equations,  respectively.  Let  us  consider  several 
influences  on  those  equations  and  therefore  on  the  vertical  profile  of  C^. 

1.  Time  changes 


During  its  early  phases,  at  least,  the  nocturnal  PBL  must  be  imsteady.  Near  sim- 
down  the  changing  surface  energy  balance  causes  the  surface  heat  flux  to  change  from 
positive  to  negative;  in  a  horizontally  homogeneous  case  the  mean  potential  temperature 
gradient  in  the  air  above  must  then  evolve  according  to 


d  de  d^w 
dtdz^  dz^  ~ 


t42j 
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In  the  initial  stages  of  evolution  dQjdz  grows  in  time,  so  that  d^6w/dz^  is  negative.  Let 
us  now  estimate  the  departures  from  a  linear  0u>  profile  this  causes. 


For  illustrative  purposes  let  the  ffw  profile  in  the  evolving  nocturnal  PBL  be 

=  Qo(l  -  z/h)  +  4SSwj-  (l  -  , 


(43) 


where  66w  is  the  midlayer  departure  from  the  linear  profile  in  steady  state.  We  scale  the 
terms  in  the  mean  temperature  evolution  equation  as 


d^Sw 

dz^ 


SSffw 

~h^' 


dt  dz 


1  A0 

T,  h 


(44) 


where  A  denotes  the  change  across  the  boundary  layer  and  is  a  time  scale.  Thus,  the 

balance  is  _ 

89w  hAQ 

Qo"  ^  87^' 

For  typiczd  values  h  =  200  m,  A0  =  5  K,  r,  =  12  h,  and  Bws  =  —0.01  m  K  s“\  (45) 
yields  ~  —1/4.  Then  to  a  good  approximation  the  Bw  profile  (43)  becomes  in 

thi«  case 

Bw  ~  0u;a(l  -  z/h)  -  Bw»^  ~  t)  =  -  z/h)^.  (46) 


The  data  on  Bw  from  the  rapidly  evolving,  early  evening  boundary  layers  in  the  Min¬ 
nesota  experiment  (Caughey  et  al.,  1979)  are  fit  well  by  (46). 

The  departure  from  linezirity  of  Bw  caused  by  time  changes  has  a  strong  effect  on 
the  Cj.  profile.  Eq.  (36)  shows  that  a  linear  Bw  profile  gives  increasing  as(l  — 
z/h)~^^^,  while  a  parabolic  Bw  profile  yields  decreasing  as  {1  —  z/h)^^^. 

Under  unsteady  but  horizontally  homogeneous  conditions  the  r-derivatives  of  mean 
horizontal  momentum  equations  are 


d'^uw  _  dV  d  dU 
dz^  dz  dtdz' 


dU  d  dV 

dz"^  dz  dt  dz 


(47) 


These  show  that  time  changes  modify  the  existing  curvature  of  the  stress  profiles. 


2.  Baroclinity 

Let  the  flow  be  steady  and  horizontally  homogeneous  except  for  mean  temperature 
gradients  in  x  and  y.  The  vertical  gradients  of  Ug  and  Vg  are  given  hydrostaticalty  by 


dz 


/ray’ 


aw 


dz 


g_  _  9  dT 

fT  dx' 


(48) 


20 


The  mean-gradient  equations  are  now 


d^wO  fT  fdVdUg  dU  \ 

dz  dz  dz  dz  )' 

showing  that  baroclinity  can  also  induce  curvature  in  the  flux  profiles.  From  (49)  the  ef¬ 
fect  on  the  stress  profiles  will  be  significant  if  the  geostrophic  wind  shear  is  an  apprecia¬ 
ble  fraction  of  the  wind  shear.  From  (48)  we  see  that  in  midlatitudes  a  horizontal  tem¬ 
perature  gradient  of  3  K  per  100  km,  not  an  unusually  large  value,  causes  a  geostrophic 
wind  shear  of  10  m  s“^  per  kilometer.  This  is  comparable  to  the  wdnd  shear,  so  we  con¬ 
clude  that  typical  magnitudes  of  baroclinity  do  have  a  substantial  effect  on  the  stress 
profile. 

Let  us  use  (50)  to  estimate  the  geostrophic  wind  shear  required  to  generate  signif¬ 
icant  departure  of  the  dw  profile  from  linearity.  Denoting  the  magnitude  of  the  mean 
wind  as  5,  we  scale  the  terms  in  (50)  as: 


UwS  fT  AS  AG 
^  g  h  h  ’ 


(51) 


Writing  (51)  as  an  expression  for  the  fractional  cheinge  in  midlayer  temperature  flux,  we 
find 


26we  _  fT  AS  AG 

dwg  4^  dWf 

For  midlatitudes,  €wg  =  —0.01  m  s~^  K,  eind  AS  AG  =  10  m^  s~^,  (52)  gives  a  70% 
change  in  the  midlayer  temperature  flux. 


3.  Terrain  slope 


Consider  a  noctmnal  PEL  that  is  homogeneous  in  plemes  parallel  to  a  land  surface 
that  slopes  at  an  angle  ^  to  the  horizontal.  Then  the  steady  mean  wind  gradient  equa¬ 
tions  in  coordinates  aligned  with  the  sloping  terrain  are  (Brest  and  Wyngaard,  1978) 


d'^uw  dV  gdQ,^,  d^vw 

dz^  ~a?' 


,9V  gd@  . 


(53) 


Here  7  is  the  angle,  measured  counterclockwise,  from  the  fall-line  vector  (the  vector  per¬ 
pendicular  to  the  contour  lines  and  pointing  down  the  slope)  to  the  i-axis.  Let  us  use 
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(53)  to  estimate  the  threshold  terrain  slope  that  gives  a  drainage  term  that  appreciably 
changes  the  stress  profile.  We  do  this  simply  by  equating  the  magnitudes  of  the  drainage 
and  Coriolis  terms.  This  yields 


fAST 
gA0  • 


(54) 


For  AS  =  10  m  s"^  and  A0  =  5  K  this  gives  d  ~  1-5  x  10"^,  a  terrain  slope  of  1.5 
parts  per  thousand.  This  is  the  slope  of  the  unusually  flat  Minnesota  site  (Caughey  et 
al.,  1979).  We  conclude  that  terrain-slope  effects  on  the  stress  profile  are  apt  to  be  sig¬ 
nificant. 


4.  Internal  gravity  waves 

De  Baas  and  Driedonks  (1985)  observed  nocturnal  PBL  structure  from  the  200-m 
tower  at  Cabauw,  The  Netherlands.  On  occasion  they  detected  Kelvin-Helmholtz  wav'es 
with  frequencies  in  the  range  3  x  10“^  Hz  to  2  x  10~^  Hz  that  persisted  as  quasi-steady 
oscillations  for  as  long  as  one  hour  without  breaking  into  turbulence.  Their  presence  was 
seen  not  only  in  the  time  series  of  fluctuations,  but  also  in  the  increase  of  the  vertical 
velocity  and  temperature  variances  with  height. 

Einaudi  and  Finnigan  (1981)  analyzed  a  nearly  monochromatic  wave  of  240  s  pe¬ 
riod  detected  in  a  weakly  stable,  nocturnal  boundary  layer  at  the  Boulder  Atmospheric 
Observatory.  They  found  that  linear  wave  theory  described  the  wave  properties  and  the 
vertical  structure  of  the  w  and  6  fluctuations  well.  In  a  second  paper  (Finnigan  emd  Ein¬ 
audi,  1981)  they  showed  that  while  the  wave-like  fluctuations  in  Reynolds  stress  were 
larger  than  the  background  turbulent  stress  levels,  they  were  nearly  in  quadrature  with 
the  gradients  of  wave  velocity,  so  they  did  not  extract  energy  from  the  wave  at  an  appre¬ 
ciable  rate;  as  a  result,  the  wave  could  persist  in  a  quasi-steady  state.  They  found  that 
the  turbulence  time  scale  was  longer  than  the  wave  period,  however,  so  the  turbulence 
was  not  in  equilibrium  with  the  large,  wave-frequency  fluctuations  in  shear  production. 

Coulter  (1990)  foimd  that  the  turbulence  levels  in  the  nocturnal  PBL  over  gently 
rolling  terrain  in  northeastern  Illinois  were  quite  different  on  consecutive  nights  with 
similar  mean  conditions.  Turbulence  levels  were  low  on  the  first  night,  but  during  the 
second  night  the  signals  showed  evidence  of  waves  of  two-minute  period  that  grew  and 
extended  beyond  the  PBL  depth  at  approximately  60-minute  intervals.  The  records 
showed  clearly  defined  active  and  inactive  periods,  with  the  differences  in  turbulence 
parameters  and  e  ranging  from  a  factor  of  2  to  4  between  them.  Coulter  inter¬ 

preted  the  two-minute  periodicities  as  Kelvin-Helmholtz  waves  made  detectable  to  the 
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acoustic  sounder  through  their  modulation  of  the  local  turbulence  structure.  He  inferred 
that  they  grew  rapidly  in  amplitude  with  time  until  they  overturned  and  broke  into  tur¬ 
bulence;  that  destroyed  the  wind  shear  and  temperature  gradient  temporarily,  but  within 
about  30  minutes  the  cycle  was  renewed. 

These  studies  indicate  that  the  nocturnal  PEL  can  have  strong  gravity-wave  activity 
that  modifies  its  turbulence  structure.  In  view  of  the  Finnigan-Einaudi  finding  that  the 
time  scales  of  the  wave  and  the  turbulence  can  be  of  the  seune  order,  we  must  conclude 
that  the  local-scaling  relations  for  structure-function  parameters  could  be  modified  as 
well.  For  this  reason,  in  his  analysis  Nieuwstadt  (1984)  excluded  cases  with  strong  grav¬ 
ity  waves  by  including  only  those  for  which  u’^  decreased  continuously  with  height.  He 
then  used  a  high-pass  filter  on  the  time  series,  as  did  Caughey  et  al.  (1979). 


C.  Implications  for  remote  sensing 

We  have  shown  that  the  sensitivity  of  structure-parameter  profiles  in  the  stable 
PEL  to  typical  magnitudes  of  flow  unsteadiness,  beiroclinity,  and  terrain  slope  precludes 
their  universality  in  the  typical  nocturnal  PEL.  This  is  consistent  with  the  findings  of 
Cuijpers  eind  Koshiek  (1988).  Their  meeisured  profiles  in  the  nocturnal  boundary 
layer  were  not  fit  well  by  the  steady  state  models  of  Nieuwstadt  (1984,  1985)  or  Sorbjan 
(1986).  The  numerical  model  of  Duynkerke  and  Driedonks  (1987),  which  uses  closures 
similar  to  those  of  Nieuwstadt  but  cdlows  for  time  evolution  of  the  boundary  layer,  gave 
the  best  results. 

Nonetheless,  the  similarity  expressions  (36)  that  relate  local  values  of  the  structure 
parameters  to  local  turbulence  parameters  are  consistent  with  observations,  LES  results, 
and  model  predictions.  This  heis  interesting  implications  for  remote  sensing.  From  the 
second  of  (36)  we  see  that  a  measurement  of  yields  the  temperature  flux  directly: 


kg 


(55) 


A  measurement  of  as  well  then  yields  the  stress  magnitude  from  the  first  of  (36): 


0.055r  {clf 

kg 

Neff  and  Coulter  (1986)  and  Gossard  (1992)  have  discussed  the  measurement  of 
structure-parameter  profiles  with  groimd-based  remote  sensors. 


(56) 
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If  measured  C|-  and  have  errors  e  (Cj-)  and  e{Cl),  then  from  (55)  and  (56) 
the  errors  in  the  inferred  fluxes  are,  to  first  order, 

e{we)  _  2e{Cl)  ejr)  1  f(Cf) 

^  2  Cl  '  T  ~  Cl  2  Cl  ' 

showing  that  the  technique  is  considerably  more  sensitive  to  errors  in  Cl  than  in  C|^. 

5.  Conclusions 

The  local-scaling  hypothesis  holds  that  over  most  of  the  nocturnal  boundary  layer 
the  turbulence  is  near  its  critical  Ricardson  number  and  its  length,  velocity,  and  tem¬ 
perature  scales  are  determined  by  the  local  stress  and  buoyancy  flux.  Data  from  sev¬ 
eral  field  experiments  and  the  Mason- Derbyshire  large-eddy  simulation  study  show  that 
the  structure-function  parameters  for  temperature  and  velocity  are  consistent  with  local 
scaling.  The  vertical  profiles  of  the  structure-function  parameters  are  apt  to  be  highly 
variable  because  of  the  sensitivity  of  the  stress  and  temperature  flux  profiles  to  unsteadi¬ 
ness  of  the  flow,  baroclinity,  and  terrain  slope.  The  local  scaling  result  does  offer  the 
prospect  that  the  stress  and  temperature  flux  profiles  can  be  inferred  from  measure¬ 
ments  of  the  structure-function  parameter  profiles,  however. 
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Figure  Captions 

1.  Experimental  data  on  C-^  in  the  stable  surface  layer,  presented  in  M-0  coordi¬ 
nates. 

2.  data  from  the  Mason-Derbyshire  large-eddy  simulation  study  of  the  nocturnal 
boundarj’  layer,  presented  in  local-scaling  coordinates. 

3.  Cl,  data  from  the  Mason-Derbyshire  large-eddy  simulation  study  of  the  nocturnal 
boundary  layer,  presented  in  local-scaling  coordinates. 

4.  A  test  of  the  local-scaling  prediction  that  the  dimensionless  Cj  approaches  a 
constant.  The  Minnesota  data  are  averages  over  a  few  measurements;  the  symbol  rep¬ 
resents  the  mean  and  the  bracket  indicates  plus  and  minus  one  standard  error. 

5.  A  test  of  the  local-scaling  prediction  that  the  dimensionless  approaches  a  con¬ 
stant. 
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